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Abstract —Directional antennas and beamforming can sig¬ 
nificantly improve point-to-point wireless links when perfectly 
aligned. In this letter we investigate the extreme opposite where 
antenna orientations and positions are chosen at random in 
the presence of Rayleigh fading. We show that while the 1- 
hop network connectivity is deteriorated, the multihop routes 
improve, especially in the dense regime. We derive closed form 
expressions for the expectation of the 1-hop and 2-hop degree 
which are verified through computer simulations. We conclude 
that node density does not greatly affect the number of hops 
required between stations whilst simple random beamforming 
schemes do, thus returning substantial network performance 
benefits dne to the existence of shorter mnlti-hop paths. 



I. Introduction 

Wireless ad hoc networks are typically equipped with 
multihop relaying and signal processing capabilities and find 
application in sensor and mobile systems such as smart grid, 
industrial and environmental monitoring, search-and-rescue 
operations etc. Commonality in many of these applications 
arises in that the number and distribution of nodes in the 
networks is often random motivating the study of random 
geometric graphs [1]. In the simplest case, these graphs consist 
of a large number of points scattered in a region of space and 
paired up whenever their mutual distance is less than some 
scalar value tq. A plethora of generalizations of this basic 
model have been put forward in an attempt to understand 
the connectivity properties of ad hoc networks and suggest 
improved network design and deployment methodologies [2]. 

In random graph models, one is typically interested in 
statistical properties such as the hop distribution portraying 
the multi-hop connectivity as popularized in the theatrical play 
“Six Degrees of Separation”. In wireless ad hoc networks, 
the hop distribution was originally investigated by Chandler 
[3] who concluded that “The node density does not greatly 
affect the number of hops required between two nodes, but 
has a much greater effect on whether a connection can be 
made at all”. Recently, a recursive formula for calculating the 
probability that a pair of nodes (i,j) separated by a distance 
Tij is connected in at least fc-hops was given in [4]. The effects 
of randomly oriented directional (or anisotropic) antennas or 
random beamforming schemes (where nodes beamform in a 
random direction) were numerically studied in [5], [6] suggest¬ 
ing that such simple schemes can lead to network performance 
gains in routing, end-to-end delay, reachability, interference 


Fig. 1. Left: Example network realization with N = 100 random directional 
nodes in a disk domain, using p = /3 = e = 1 and ri = A. The blue arrows show 
the directional links available that can be utilised for routing and scheduling 
transmissions. Right: Cardioid radiation patterns for e = 1/2,1. 


tolerance, and capacity. Recent analytical work has shown 
that the reported benefits are actually highly dependent on the 
prevailing channel conditions described by path loss exponent 
T] experienced in the propagation medium [7], [8]. 

In this letter we analytically quantify these benefits for 
the first time and obtain closed form expressions for the 
multihop reachability of random networks with randomly 
oriented antennas in the presence of Rayleigh fading. We 
therefore show when and how schemes like random beam¬ 
forming provide progressively better multihop coverage and 
reachability (measured by the number of accessible nodes in 
at least k hops). This follows from the more likely existence of 
shorter multi-hop paths between stations in the dense network 
regime with randomly oriented directional antennas. Implicit 
for this performance improvement is perfect interference man¬ 
agement and a good MAC. We derive communication theoretic 
lower bounds on the multihop connectivity of networks with 
randomly oriented directional antennas, useful when designing 
wireless sensor networks (WSNs), for instance, for choosing 
the right density of directional nodes to be deployed in order 
to meet certain multihop requirements or performance goals. 
Significantly, we also confirm Chandler’s conclusion [3] by 
showing that the typical hop distance h ^ where the 

coefficient c strongly depends on the path loss exponent rj and 
antenna directivity. The exact form of c remains an interesting 
open question. 
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II. Network Definitions and System Model 

Consider a wireless network with identical nodes located 
randomly on the plane according to a spatial Poisson 
point process (PPP) of intensity p. We denote the location 
of node i by the vector = {ri,9i) in polar coordinates. 
Furthermore, consider the case where each node is equipped 
with a directional antenna pointing in a direction -di chosen 
randomly in [0, 27r]. Such a configuration is commonly found 
in WSN applications where sensors or smart meters form a 
random mesh topology (see for example Fig. 1). 

Assuming negligible inter-node interference (e.g. perfect 
CDMA/TDMA), we define the connection probability between 
transmitting node i and receiving node j through the relation 
Hij = P(SNR • |/ip > p), where SNR denotes the long-term 
average received signal-to-noise ratio and h is the channel 
transfer coefficient for single input single output (SISO) an¬ 
tenna systems. Assuming identical lossless antennas we have 
from the Friis transmission formula that SNR c>c GijGjir~p, 
where Vij = Ir—rj| is the Euclidean distance between the two 
nodes, p is the path loss exponent (typically p > 2 ), and 
is the gain of the antenna at node i observed in the direction of 
node j. Notice that Hij = Hji since the channel is assumed 
reciprocal. Isotropic antennas have a constant gain G = 1, 
while anisotropic ones are functions of the polar angle 0, 
appropriately normalized by the condition G(0)d0 — 27r. 
Assume that the antenna gain is given by a cardiod function 

Gij(9ij) = 1 -f e cos 9ij, ( 1 ) 

where e G [ 0 , 1 ] measures the extent of deformation from the 
isotropic case (e = 0 ), and 9ij is the direction of receiving 
node j relative to the antenna orientation of node i (see Fig. 
1). The cardiod function (1) is representative of wide-angle 
unidirectional radiation patterns often found in microstrip 
(patch) antennas. We expect however that the results presented 
herein qualitatively apply to other directional gain functions 
as well as multi-directional e.g. G = 1 -I- e cos n9 for n > 1 
as was the case in previous studies [ 8 ]. 

The framework that follows is statistical and assumes a 
Rayleigh fading channel such that |/ip ~ exp(l) and 

= P (|/ip > -p/SNR) = (2) 

where /3 depends on the transmission wavelength, signal 
power, the threshold p etc., and defines an effective communi¬ 
cation range rg = [GijGji/. In simulations we set /3 = 
1. Connection probability functions like (2) have been used 
extensively in the literature however one should note that there 
are many other possible functions (e.g. including shadowing 
or a different path loss function) and it is not clear which one 
models real networks more accurately. Equations (1) and (2) 
are therefore examples used for the sake of juxtaposition with 
some of our results and approach being more general. 

III. K-hop Connectivity 
A. Mean 1-hop degree 

It is often desirable in WSNs that broadcasts are heard by 
as many neighbouring nodes as possible. A simple way of 


measuring this is the mean degree. Denoting the 1-hop degree 
of transmitting node i by d^, the 1 -hop mean degree pi is 
given by the arithmetic mean of Equivalently, since nodes 
are Poisson (i.e. uniformly and independently) distributed, 
we take the spatial and orientation average of equation ( 2 ) 
multiplied by the density and normalised by 2 ?! and calculate 


/•27r /•2'7r ^oo 

^ Jo Jo 


pr(2/77) 

27r?7/3^/’i 



(3) 


In the first line of (3) we have used the homogeneity of 
the spatial PPP on and set node i at the coordinate 
origin, oriented along the x-axis. Therefore, in (3) the first 
two integrals {drjd9j) average over all possible positions of 
Tj G K^, and the third integral (di?j) averages over all possible 
orientations dj G [0,27r] of antenna j. In the last line of (3) 
we have used the periodicity property of the cardiod function 
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where 2 F 1 is the Gauss hypergeometric function. Eor any 77 > 
2, Pi is monotonically decreasing with e. This implies that 
the more directional the antenna gain is, the less nodes can 
be reached in a single hop. Note that since we have assumed 
identical antennas for both transmitter and receiver chains (i.e. 
Hij = Hji), the in- and out-degree of network nodes are equal. 


B. Mean 2-hop degree 

The 1-hop mean degree pi is a good measure of the local 
connectivity of a network. Often however, and particularly in 
WSN MAC design [9], it is desirable to maximize the number 
of 2-hop neighbours that a broadcast can reach. Let the 2-hop 
degree of transmitting node i be denoted by d\ counting the 
number of 2 -hop neighbours of node i which are not 1 -hop 

neighbours. To make progress, define the probability of nodes 

( 2 ) 

i and j having at least one common 1 -hop neighbour as H^ 
which can be expressed as the complement of not having any 

=1-n (5) 

where nodes k ^ i,j act as a relays. Using (5) and assuming 
that |ri| = = 0 we may express the 2 -hop mean degree as 

^27r n‘2'K pOO 

^d){HS'')rjdrjd9jd'0j, (6) 

(7T^^^)=1 - exp J HikHkjrkdrkddkddk^ (7) 

It is therefore immediately clear that p 2 is not a linear function 
of the density p (c.f (3)). Eor example, for a simplified 
unit disk model (i.e. ignoring fading) it can be shown (see 
Appendix I.) that p 2 ~ ^P'^t'o ~ 27r(2ro)^^^r(2/3)(p/3)^/^. 
Does the mean 2-hop degree p 2 increase or decrease with 
antenna directionality e? This cannot be answered by simply 
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looking at equation ( 6 ) and therefore we postpone further 
discussion for the next section where we numerically evaluate 
the integrals and compare with computer simulations. 


C. Mean k-hop degree 

Generalizing ( 6 ) to the case of the fc-hop mean degree is not 
straight forward. Let the fc-hop degree of transmitting node i 
be denoted by d\ counting the number of fc-hop neighbours 

of node i which are not {k — l)-hop. The fc-hop mean degree 

ik) 

fik is therefore given by the arithmetic mean of . Define 
the probability that nodes i and j have at least one common 
TO-hop neighbour, which with some care can be expressed as 


«r’=1 - n 


l-H. 


(m-l) jj(l) 


ik 


kj 


m—z 

n( 




( 8 ) 


where we have used = Hij for ease of notation. Equation 
( 8 ) is a nested equation calculating the complement probability 
of two nodes having no common m-hop neighbour, which is 
itself given by a similar expression. The fc-hop mean degree is 
therefore given by ^ times the expected probabilities of two 
nodes being fc-hop neighbours but not m-hop, for m < fc — 1 


^ fc —1 

n ■ -drNMN (9) 

in the limit of the number of nodes N —>■ oo and where we 
have dehned Ad = x [0, 27r]. Equation (9) is not particularly 
insightful and cannot be easily computed, even numerically, or 
even after some approximations [4]. In analogy to the result in 
Appendix I. one may expect that pk = p(2fc—l) 7 rrg— 

We will conhrm this in Sec. V for the case of fc = 3, however 
are unsure of the correction term for higher values of fc. 


IV. The Typical Hop Distance 


Given a large but hnite network of N nodes, pick two 
nodes at random. What is the expected hop distance h between 
these two nodes? Numerical studies have shown that random 
beamforming schemes can reduce this number in mesh net¬ 
works, thus resulting in faster message dissemination while 
also reducing the total signal processing done by relays [5]. 
To provide a formal definition of the typical hop distance 
we begin from the identity stating that the sum of all fc-hop 
neighbours of any node i equals to — 1. Symbolically this 
is: d[°°^ + d!'^'^ = N — where is the number 

of nodes which cannot be reached by i in any number of 
hops. Note that 1 -f connected 

component that node i belongs to. In fact, averaging this over 
all nodes gives the typical cluster size in the network. Clearly 
in a fully connected network d^^'^ = 0, Vi. Averaging both 
sides of the identity over i and changing the summation order 

N / N-1 \ N-1 

i=l \ k=l ) k=\ 

indicating that p,k/{N — 1 ) is the probability mass function 
of the hop distribution for large random networks. Signih- 
cantly, the typical hop distance h between a random pair of 



Fig. 2. Left: Simulation results for the fc-hop mean degree as a function 
of the density p for the isotropic (e = 0) and anisotropic (e = 1) case using 
77 = 3. For fc = 1, the two black straight lines are given by (3). For fc = 2 
the black curves are obtained by numerically integrating (6), while for fc = 3, 
the black curves are fitted to a — bp^^^ + cp for a,b,c > 0. The markers 
are obtained from computer simulations in a large network similar to that in 
Fig. 1. Middle and Right: Parameter space (p,7?) showing isotropic (red) or 
anisotropic (blue) radiation patterns deliver higher pi and p 2 respectively. 


nodes is given by the mean of this hop distribution h = 
]vE ^dk- We have therefore established a link between 

the typical hop distance h and the fc-hop mean degree (9). 
Notice that if the network is not fully connected (i.e. 0) 

then h diverges. 


V. Numerical Analysis 

In the simulations that follow, we independently and uni¬ 
formly generate the spatial and orientation coordinates of 
N = YpV\ nodes in a circular domain V of area V = ttR^ and 
radius i? = 10 using /3 = 1. Whenever a randomly generated 
number C, G [0,1] < Hij, nodes i and j are paired up (see 
Eig. 1). This guarantees that the network links are statistically 
independent. The resulting undirected graph connections are 
stored in a symmetric adjacency matrix from which we extract 
useful network observables such as the fc-hop degree of each 
node. In order to keep our results free of boundary effects, we 
use values drawn from nodes situated away from the border 
of V as done in [10]. Einally, we improve our statistics by 
repeating the above process in a Monte Carlo fashion. 

In the left panel of Eig. 2 we plot the simulation results 
for pk for fc = 1, 2, 3, as functions of the node density p for 
f] = 3. Pi increases linearly and the isotropic case is slightly 
higher (10%) than the anisotropic one. The two straight lines 
are plotted using (3). p 2 and p^ increase super-linearly and the 
anisotropic case is higher by 20% and 35% than the isotropic 
one at p = 4 respectively. Eor p 2 , the two curves are obtained 
by numerically integrating (6), and for p^ they are fitted to 
a — bp^^^ + cp for constants a, 6, c > 0 

The middle and right panels of Eig. 2 shows in different 
colors the gain (isotropic or anisotropic) which maximizes p^ 
in the parameter space of (p, p) for fc = 1,2. Looking at the 
middle panel of Eig. 2, it is clear that for p > 2 networks with 
isotropic antennas will typically have a higher 1-hop mean 
degree than with anisotropic radiation patterns. Signihcantly, 
this characterization is independent of the density p. Eor fc = 2 
however, anisotropic patterns are superior and occupy more 
and more space at higher node densities. We therefore reach 
the following non-trivial conclusion: multihop accessibility is 
intensified for directional antennas, especially at high node 
densities. We expect that the more directional G is, the better 
its multi-hop accessibility will be in the dense regime. The 
mechanism by which this occurs is simple: 1) Fixed p: as 
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Fig. 3. Left: Comparison of the hop distribution between isotropic (blue 
bars) and anisotropic (purple bars) networks using p = r] = 3. The y-axis 
therefore gives the probability that two randomly selected nodes are fc-hops 
apart. Right: Plot showing the typical hop distance /i. as a function of the node 
density p for different values of (e,??) with fitted black curves for ~ . 


node density p increases, the next hop node is more likely 
to be oriented in a better direction (than from the forwarding 
node), resulting in better performance (in terms of reachability) 
of networks with anisotropic antennas. 2) Fixed p: high path 
loss exponent p degrades the performance of both radiation 
patterns, but does so more for the anisotropic one. This is 
because concentrated beams are more likely to be obstructed 
in cluttered environments resulting in better performance of 
networks with isotropic antennas. 

Fig. 3a. shows the hop distribution obtained from a col¬ 
lection of networks generated using p = rj = 3 thus giving 
further insight into the comparison of reachability and network 
connectivity for isotropic and anisotropic antennas. Recall 
that each bar in the histogram corresponds to the mean k- 
hop degree p,^. As expected, pi is slightly larger for the 
isotropic case (blue) since p > 2. For A: > 1 the anisotropic 
distribution is skewed to the left indicative of the reachability 
benefits of networks with anisotropic antennas. Physically 
this means that the resulting network has shorter paths (in 
terms of hops) between pairs of nodes; a clear indicator 
of potential gains in end-to-end delay subject to efficiently 
designed MAC and routing protocols [9]. Fig. 3b. shows plots 
of the typical hop distance h{p) for different values of (e, p). 
h rapidly increases to a maximum at which point the network 
is fully connected, and then slowly decays like ~ thus 

conhrming Chandler’s conclusion [3]. It is clear however that 
the typical hop distance h is about 15% less in mesh networks 
with randomly oriented directional antennas. 

VI. Conclusions and Discussion 

In this letter we have studied the multihop connectivity of 
networks with directional antennas. Specifically we have fo¬ 
cused on random beamforming strategies where nodes choose 
the direction of radiation at random; a practical solution when 
no a priori node location information is available or in random 
uncontrolled (e.g. air-dropped) WSN deployments. We have 
developed a general mathematical framework and investigated 
the interplay between node density p, antenna directivity G, 
and the pathloss exponent rj. To this end, we have shown 
that random beamforming strategies can substantially improve 
multihop coverage, reachability, and also reduce the typical 
hop distance between nodes, thus improving end-to-end delay 
and signalling overheads. Signihcantly, we have obtained 
closed form expressions for the expected one- and two- hop 
mean degree thus enabling further design optimization of 


large wireless networks with randomly oriented directional 
antennas. Furthermore, knowledge of such statistics can be 
used to make MAC and routing protocols smarter. In the 
future, it would certainly be interesting to conhrm our hndings 
using packet-level network simulators, e.g. ns3, or even in real 
experimental test-bed deployments. It would further be inter¬ 
esting to extend this theoretical work on randomly oriented 
directional antennas to interference limited wireless networks 
by means of stochastic geometry methods. 

Appendix I 

2-hop mean degree for the hard disk model 
Consider the hard disk connection morel given by Hij = 1 
if Tij < To and 0 otherwise and substitute it into ( 6 ) to get 

P 2 = 27rp / Tij fl - drij (I.ll) 

Jrg ^ ^ 

where A(ry) = 2 rQ arccos ^ — rT is the in¬ 

tersection area of two overlapping circles of radius tq, a 
distance Vij apart. As A{rij) appears in the exponential of 
(I.ll), the dominant contribution to p 2 at high node densities 
will be due to a situation when there is little to no overlap 
between the communication regions of nodes i and j i.e. when 
A{rij) 1 and a common 1-hop neighbour is unlikely. We 
therefore expand A{rij) near « 2 ro and obtain A(rij) « 
-^^^(2ro — Substituting this into (I.ll) we get 

P2 « 27rp / rij ^1 - dr^ 

Jtq ^ ^ 

= 3pTtrl - 2^(2ro)2/3r (2/3) (p/3)'/" + 0(p-'/"), 
since A{rij > 2ro) = 0. In the 3D network case we get 
P2 ^ 28/37rpr^ — 87rrg/^-y/2p. 
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